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ABSTRACT 

We consider the effects of strong critical layer nonlinearity on the 
spatial evolution of an initially linear "acoustic mode" instability wave on a 
hypersonic flat plate boundary layer. Our analysis shows that nonlinearity, 
which Is Initially confined to a thin critical layer, first becomes Important 
when the amplitude of the pressure fluctuations become 0(1 /M 4 In M 2 ), where 
M Is the free stream Mach number. The flow outside the critical layer is 
still determined by linear dynamics and therefore takes the form of a linear 
Instability wave - but with Its amplitude completely determined by the flow 
within the critical layer. The latter flow Is determined by a coupled set of 
nonlinear equations, which we had to solve numerically. 

1. INTRODUCTION 

Laminar boundary layer Instabilities are predominantly invlscid at 
sufficiently high Mach numbers (Mack, 1984, 1987) with the so-called vorticlty 
modes exhibiting the most rapid growth at very large Mach numbers. The 
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so-called acoustic modes exhibit the most rapid growth at more moderate Mach 
numbers - with the two-dimensional disturbances growing more rapidly than the 
corresponding oblique waves. 

Mack (1984, 1987) computed the relevant numerical solutions to Rayleigh's 
equation and Cowley and Hall (1988) worked out the corresponding asymptotic 
solution for the hypersonic limit where the free stream Mach number M -» ®, 
Their results, as well as those of Mack (1984, 1987), suggest that, while the 
Instability wave number becomes small, the instability wave growth rate becomes 
even smaller as M -* « which means that there Is a well-defined critical layer 
at every unstable frequency in this limit. Nonlinear effects will then balance 
the resulting singularity at sufficiently large Reynolds numbers and the 
present work Is concerned with extending the Cowley-Hall (1988) analysis into 
this nonlinear regime. 

Boundary-layer-transition experiments often Involve spatially growing 
Instability waves generated by relatively two-dimensional, single-frequency 
excitation devices such as vibrating ribbons or acoustic speakers. As the 
Instability wave propagates downstream, its amplitude continues to Increase 
until nonlinear effects come Into play - provided, of course, that the Initial 
amplitude Is sufficiently large and/or the mean flow divergence Is sufficiently 
small (i.e., the Reynolds number Is sufficiently large). 

Our previous remarks suggest that the nonlinearity first comes Into play 
locally, I.e., within a critical layer. The solution to the resulting 
nonlinear critical layer problem reduces to the linear small growth rate 
hypersonic instability wave solution far upstream, and Is therefore the natural 
downstream continuation of this upstream linear solution. The flow outside the 
critical layer continues to be linear but the corresponding Instability wave 
amplitude Is completely determined by the nonlinear dynamics with the critical 
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layer. These considerations fix the relative scaling between instability wave 
amplitude and Mach number and thereby show that the nonlinearity first becomes 
important when the pressure fluctuation amplitude becomes 0(1 /M^ In M^) In the 
main boundary layer. 

The critical layer nonlinearity is strong in the sense that it enters 
through a coefficient in the lowest order equation (as in the Goldstein, 

Durbin, and Leib, 1987 incompressible boundary layer analysis) rather than 
through an inhomogeneous higher order term (as in the Goldstein and Leib, 1989 
compressible shear layer analysis). However, it differs from the former 
analysis (but is similar to the latter) in that coupled vorticity and energy 
equations have to be solved simultaneously. This has been done numerically. 

While nonlinear effects always decrease the instability wave growth in 
the Goldstein et al . (1987) analysis, compressibility effects produce a 
dramatic Increase in growth in the present analysis - leading to a possible 
singularity at a finite downstream position (as found by Goldstein and Leib 
(1989)). Finally it is worth noting that there Is an unusual nonlinear 
modulation of the basic critical layer streamline pattern which does not 
usually occur in nonlinear critical layer analyses. 

The overall plan of the paper is as follows. The problem is formulated 
in section 2, where we show how the nonlinear flow gradually evolves from the 
strictly linear hypersonic solution. The flow outside the critical layer is a 
linear invlscid perturbation about a hypersonic (i.e., M >> 1) boundary layer 
flow, and is found by extending the asymptotic analysis of Cowley and Hall 
(1988) into the nonlinear regime. This flow has a double layer structure and 
the solution for the main boundary layer region is worked out in section 3. 

The critical layer is contained in the outer region and in section 4 we obtain 
the outer region solution that applies outside the critical layer. This is 
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then used to formulate the relevant critical layer problem in section 5. The 
resulting critical layer vorticity and energy equations are different from the 
usual nonlinear nonequilibrium transport equations but can be transformed into 
the latter by a suitable change of Independent variables. The numerical 
solution of these coupled equations is described in section 6. Finally, 
numerical results are discussed in section 7. 


2. FORMULATION 


The free stream flow parameters are used as reference quantities and are 
generally denoted by the subscript ®. The reference length, say A is taken 
to be some suitable boundary layer thickness (e.g., momentum thickness). Then 
the steady flow is characterized by the Mach number 

M = Uc/Coo (2.1) 

and Reynolds number 


Re E UcoA/v®, (2.2) 

where 

Co - (y^T*,)!^ (2 .3) 

Is the speed of sound in the free stream, u is the kinematic viscosity, y is 
the isentropic exponent of the gas, and & is the gas constant. 

We suppose that the flow is two dimensional and that Re is large enough 
that the unsteady motion is essentially inviscid and unaffected by mean 
boundary layer growth over the region in which nonlinear interaction takes 
place. We can then suppose that the mean pressure is constant and that the 
mean flow velocity U(y) depends only on the transverse coordinate y to the 
required order of approximation. We assume that the wall is Insulated and, 
for simplicity, restrict our attention to an ideal gas which satisfies 
Southerland's viscosity law and has Prandtl number unity. In this case the 
local mean density R and mean temperature T will be related by 
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RT = 1 , 

and the mean velocity and temperature will be given by 

U = h 1 <n> 


and 


(2.4) 

(2.5) 


T = 1 + M 2 (l - U 2 ), (2.6) 

respectively, where n is the Dorodni tsy-Howarth variable (Stewartson, 1964) 
defined by 


n 


r y si 

J o T ' 


(2.7) 


the prime denotes differentiation with respect to n» and h is the Blasius 
function, i.e., it satisfies 


hh' 1 + 2h" ' =0. 

It follows that (Schl 1 chting , 1960, pp. 117 and 118) 


( 2 . 8 ) 


U ■ 1 — — e‘ ; ' 2/4 


'-k* <0 


bV ^' 2 


(2.9) 


as n - 00 , where 


n = n - 3, (2.10) 

3 = 1 .73, and b * 0.462. 

As Indicated in section 1, we suppose that the unsteady motion starts as 
a linear Invlscid Instability wave (which is governed by Rayleigh's equation) 
in the upstream region where x - -=°. We also suppose that the linear mode is 
of the acoustic type (Mack, 1984, 1987) and that 

a E ^ << 1 . (2.11) 
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The relevant asymptotic solution to Rayleigh's equation was worked out to 
lowest approximation by Cowley and Hall (1988). They show that the scaled 
complex wave number a is 0(a), that the phase speed c behaves like 

c = 1 - oc ~ 1 as a - 0, (2.12) 

where c = 0(1), and that the flow develops a double (actually triple if the 
free stream Is included) layered structure, with a relatively thin outer 
region, where 


In Y = 6^26 - n) 

is order unity provided the constant 6 satisfies 

.2 


r s 

26 


a. 


(2.13) 


(2.14) 


Cowley and Hall (1989) extended their analysis to higher order In a to show 
that the instability wave growth rate (or imaginary part of the complex wave 


number) is small relative to i?ea, or more precisely that 



(2.15) 


The linear instability wave will then have a distinct critical layer, 
which, in view of (2.9), (2.12), and (2.13), must lie in the outer region where 


Y = 0(1), and 


U 


obY 1 



In Y 


I In 2 ? 


: ) 


0(6 3 ) 



(2.16) 


We expect the analysis to become nonlinear in this layer because the flow 
is assumed to be inviscid and the linear instability wave amplitude Increases 
in the downstream direction. The motion Is predominantly linear outside the 
critical layer and we expect the solution there to expand like 

2 

U = U(n) + eUj + e lij + . . . (2.17) 
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V 


( 2 . 18 ) 


2 

= eVj + e ^2 + • • • 

9 = T(n) + ex-j + + . . . (2.19) 

p 1/y = 1 + eir ] + e 2 ir 2 + . . . , (2.20) 

where {u,v} are the velocity components in the x and y directions, 9 is 
the temperature, p is the pressure, and e Is a characteristic amplitude of 
the instability wave in the streamwise region where nonlinearity first becomes 
important. It will be specified more precisely below. Then since 
{ui ,vi , ti , pi } are determined by linear dynamics, it follows from Goldstein 
(1984) that the expansion coefficients in (2.17) to (2.20) (which depend on x, 
n, t, and a) are determined by 


Ltr ] = 0 

(2.21) 

Lir - 2_ (M. + 2U1 93 

Lir 2 " Dt V9x + 9y/ T 3x 

(2.22) 

D U' T 8ir l 

Dt U 1 + T V 1 = “ ctT 9 x 

(2.23) 

D U ' 

Dt <U 2 + U l Tr l > + T~ (V 2 + V l ff l ) = _aT IF “ f 

(2.24) 

Dv , Sir, 

dt - -° T W 

(2.25) 

D 8ir 2 

Dt (V 2 + V l ir l ) = ~ oT IF “ 9 

(2.26) 

D 

5t T i + rv i - ( ^ 1 )t dT 

(2.27) 


and a similar equation for T 2 , where the prime still denotes differentiation 
with respect to n. 

L E 5t (^I - <’ 7 • T? ) * 2oU ' a T5? (228) 
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(2.29) 


and 



(2.30) 

(2.31 ) 


9 E §7 Vl * §7 V 1 + <”1 ly V (2.32) 

It now follows from (2.12), (2.15), and the fact that a = 0(a) that the 
solution to (2.21), (2.23), (2.25), and (2.27) that matches onto the upstream 
linear solution must be of the form 


u 1 - j ^f 1 (n,x 1 )A + (x 1 )e 1X , 


(2.33) 


= -a &1<Mn.x,)A + e iX , 


(2.34) 


■ ^2 ^®®iA^e^ X , 

a 


(2.35) 


where we have put 


TTj =. ^en 1 A + e lX , 


, A 

1 = S’ 


(2.36) 


(2.37) 


X = aa£x - (1 - aC)tj, 

A"*’( x i ) , a (a), and c(o) are real quantities and, to the required level of 
approximation, n-j satisfies Rayleigh's equation 


(2.38) 


= o, 


(2.39) 


where 
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E (U - c> 2 i X - (na) 2 lTT - - (U - c) 2 l for n * 1, 2, . . . 

an (u _ c ) 2 dr > La J 

(2.40) 

are the linear Rayleigh operators, the complex wave number and phase speed a 
and c, respectively, are given by 


a = aa + 


a? All 
6 lA*’ 


(2.41) 


2 4 +' 

_ i - a A 

C aC ” f ’ 
Sa 1A T 


(2.42) 


and the prime now denotes differentiation with respect to X] . We note In 
passing that a has an expansion of the form 


_ _ _ 2 - 
a = + a<X 2 + a + . . . , 


(2.43) 

where each of the coefficients has Its own series expansion* In , In 5, 
etc. (see (2.14)) and similarly for c, l.e.. 


-(0) 1 -d) 

a = a + 7 a + . 
n n 5 n 


(2.44) 


- -( 0 ) 1 -( 1 ) 

C = C + 7 C„ + . 

n n 5 n 


n ] must satisfy the boundary conditions 

dlL 


(2.45) 


1 


dn 


= 0 at n = 0, 


and 


n 1 - 0 as 


(2.46) 


(2.47) 


and, In order to match with the linear solution far upstream, we must require 
that 


*It might be helpful to think of a and 6 as Independent expansion 
coefficients at this point. 
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A 


t 



as x, - -®, 


(2.48) 


where k Is the scaled growth rate of the upstream linear Instability wave 
and a Is a complex constant. Finally, the remaining functions of n. I.e., 
*1. 4»i, and 0] can readily be found from (2.23), (2.25), and (2.27) once 
II] is known. 

In the following two sections we derive dispersion relations for the 
Instability wave amplitude outside the critical layer. 

3. LINEAR SOLUTION IN MAIN BOUNDARY LAYER 
First suppose that n = 0(1). It is easy to see from the Cowley-Hall 
(1988) analysis that II] should expand like 


n i = a ( P 0 + ctP 1 + + ' • •)’ 

where each of the expansion coefficients In turn has the expansion 

p p <0) A 1 D (l) 1 d (2) _ , „ 

P n " P n * 6 p n * ~2 p „ * ■ ■ ■ for n . 1. 2, . . . 

6 


(3.1) 


(3.2) 


Substituting these Into (2.39) and equating coefficients of like powers of 
a and S, we find that 


9t 



o, 


(3.3) 



where 


(0) n (l) „ (1 - U) 2 

1 k 2 " g 2 P 


(3.5) 
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(3.6) 


Q 2 = Q 2r ‘ 1q 


it • 


Q 2r q are real, 


q = (y - 1X1 + U> a] [^-f 1 (1 + U) + u]p 2 + 
and we have put 


c^Cl - U) 


u- dP 0 

4 dn ' 


(3.7) 


g\ 0) = (1 “ U)2 h 1 


,2 dn 


(U - 1) 

- _ u 2 >[*-=-^ (1 - U 2 ) - 


(U - 1)‘ 


(3.8) 


Equation (3.3) was solved numerically by Cowley and Hall (1988). They point 
out that 


P - 5- e-' 2 ' 2 
r 0 ~3 e 
n 


as 


(3.9) 


where D is an, as yet, undetermined constant. Equations (3.4) and (3.5) can 
now be solved by variation of parameters to obtain 

' ( N2 f n 

P 1 = P 0 V/ Vn>dn dr i’ (3.10-a) 


P (1> = P 
*2 K 0 


1 - U 


/N ^ 


0 2 (n)dn dp. 


It follows that 


p (0) _ 
r l D 


r 00 ( n u 2 p“ 

0, d n ; P 2 ~ g- J Q 2 dn, as n - «. (3.11-a,b) 


(3.10-b) 
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4. LINEAR SOLUTION IN THE EDGE LAYER 

Cowley and Hall (1988) point out that the expansion (3.1) breaks down at 
large disturbances from the wall and then proceed to construct a new "outer" 
solution for the region where Y = 0(1) (see (2.13)). The expansion In this 
region must be of the form 


n i = 1 + {t)P ] + (f) P„ + . 


(4.1) 


where again each expansion coefficient has Its own series expansion 


n n n S n 


(4.2) 


in terms of S - ^, In 6, etc. Substituting this, together with the new variable 
(2.13), Into (2.39) and equating coefficients of (a/S), we find upon 
integration that 


p(-l) „ |<-1> 

' 1 t l 



b' 0) (^ 


- 2c 1 bY + 



?( 0 ) 
+ E 1 


(4.3) 





i_ai: 

«1 1 A f 







J 0 



. - c? 2 - 

In Y dY + j 1 In Y - 



X 


( b? - c i) : 





(4.4) 


where we note that the lowest order solution < 1 . e . . Pj ) Is given In Cowley and 
Hall (1988), , E«, and 4» ± are constants of Integration (the latter of 
which can be different depending on whether Y5;0), and 


12 


(4.5) 


T c = 1 + (y - Dc 1 

Is the mean temperature at the critical level where 

bY = c r (4.6) 

This solution does not satisfy appropriate free stream boundary conditions and 
it Is necessary to Introduce another outer expansion for the region where the 
variable 

n = an (4.7) 

Is order one. The solution In this region Is 


n 1 = (1 + aB.| )e 


~ -a /k 


2 ^ 
ac 1 n 


(4.8) 


to the required order of accuracy and matching with (4.1) shows that 

5 ( 0 ) _ ^ 

B 1 " - 2 ’ 


1 


and 


(4.9) 


5(0) 


“lf T c " " 5 1 (y " 1>2 ] [♦" + 1n l C l l 


]) 


2_ f-(l> 1 A f,> \ 

+ Ci \ 2 + a] ia7 


(4.10) 


Finally, matching (3.1) and (4.1) and using (3.9) and (3.11) shows that 

b 2 «, 


:2 * 


(4.11) 


c 2 - 


1 = r 1 [ Q, d n , 

nr i a 1 


a i °o 


(4.12) 


Dm 4 > + ■ 0 , 


(4.13) 


and 
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(4.14) 


“i'i[' - - 1 ) 2 ]A*+" * 2 Jr ■ 5 i Jr J 0 q dn - 

Substituting (2.13), (2.34), (2.36), (2.41), (2.42), (4.3), and (4.4) into 
(2.25) and equating coefficients of like powers of a, we find 




(C - 

bY) 

* f°\ 1 

fj2 1 

aH 

( s i. 

) 2 

W 5! 

T c + _ 2 _ 
v a l C] 

ia V 


(4.15) 


Substituting this together with (2.33) and the previous equations into (2.23) 
and using (2.9) 



And finally, proceeding similarly with (2.35) and (2.27), we obtain 

Vy - 1) „ 


a, - - 


r(ac,) 2 


bY + 


(4.17) 


5. NONLINEAR EDGE LAYER TERMS 


The lowest approximation (in terms of a) of the 0(e 2 ) terms in the 
expansions (2.17) to (2.20) has to be determined before the solution within 
the critical layer can be found. However, it is only necessary to consider 
the edge layer solution for this purpose. It follows from (2.33) to (2.36) 
that this solution must be of the form 


*2 



Ren 


(2) 2iX 
2 e 


(5.1) 
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V 2 + V 1 ff l = ~^ a ^2^ 


♦ /W 2 ) e 2ix 


ii „ „ 5 *(0) /) ^<2) 2iX 

U 2 + U l Tr 1 “ a” *2 + KB ^2 e 


where the, as yet, unknown coefficients are functions of n. X] , and a. 
Substituting (5.1) into (2.22) and using (2.31), (2.28) to (2.36), and (2.40), 
we find 

S'.n' 2 ’ - - - wlaTF . (U - c> 2 *j- f 2 J } , (5.4: 

22 5 \_ dn L<u - c) 2 J J 

where 




Substituting (2.6), (2.13), (2.16), (2.40) to (2.43), and (4.1), (4.3), and 
(4.15) to (4.17) into this result we find that 


Yc, - bY 


f 1- 
J (~- 


dY (c, - bY) 2 dY acj T 


1 bY u h - bY ' 

^ r bY + — — 


bY(y - l)(c 1 - bY 


to lowest approximation in o when Y = 0(1). It follows that 


r xl (c, - bY> t2 


c, - bY 

». I j 


as Y - jp. 


(5.8) 


15 



It therefore follows from (2.13), (2.16), (2.26), (5.2), and (5.6) that 


<I> (2) * i- 
v 2 2 


i+2 


a" 4T c c^a 2 


[' ♦ °( c "i - "?)]• 


as 


11 

b ’ 


(5.9) 


and similarly from (2.24), (5.3), and (5.5) that 

, t ,( 2 ) 6 


2 [ J ! + 0 ( b7 - f ,)] 


as 


Y - — 

b ’ 


(5.10) 


where Is an 0(1) constant. 

6. THE CRITICAL LAYER 

Equation (4.16) shows that the edge layer solution becomes singular In the 
critical layer where Y = b/^ . The governing equations therefore have to be 
rescaled to obtain a bounded result In this region. The thickness of the 
linear smal 1 -growth-rate critical layer Is of the order of that growth rate 
divided by the mean velocity gradient times the real part of the wave number, 

1 .e. , 



It therefore follows from (2.10) and (2.13) that the appropriate transverse 
coordinate in this region is 

( 6 . 1 ) 

Equations (2.16) to (2.20), (2.33) to (2.36), (4.1), (4.15), (4.16), (5.1) 
to (5.3), and (5.8) to (5.10) suggest that the flow in this region should 
expand like 
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u = 1 - o( b 


, h c 2 b7 (° C A c6 /W e ,x 

' b ' 5 i 25 ) ^^7 


+ eu 


*3 2 
S e - 

* — u 2 * ■ 


( 6 . 2 ) 


i ) Re 

a 


bY 

ii 


?)■’ -?- r 

i y c i a i 


IX S 3 e 2 1 


4 5"?3 t 

a 2aC.,T 
I C 


&1A +2 e 21 * 


+ CVj + . 


e = T c + f-(y - l)bY + + . 


(6.3) 

(6.4) 


and 


p 1/x - 1 * c &aV X + f-i, ♦ . . . , 

where the passive terms Involving ln(c/6) have been Incorporated Into u ] , 
etc., and 

2 


(6.5) 


v = v — . 

\Oj 


( 6 . 6 ) 


Then the critical layer solution will match with the "outer" edge layer 
solution If we require that 


- - ^I> /? eA V x , 


T cVl 


(6.7) 


au 


^ - -tSf 

3Y S 1 (T c e,) z 


(6.8) 


as 


Y - ±<*>, 


and 
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(6.9) 


Au 


1 


HmfujCY) - (-Y)j = 

Y-oo 


1 - 


.2, 

•1 


C7(y - 


1)‘ 


-x.tJX 


T c C l 


y?e(<|> + - <J> )A T e 


+ higher harmonics. 


The expansion coefficients Op v ] , x ] , etc. are functions of X, Y, and 
x ] only. They are determined by the Invlscld vortlclty, energy and continuity 
equations, which can be written as 



and 



<Y 



T ”y 



where we have put 


and 


D = a(u 





(6.1 l-a,b) 


( 6 . 12 ) 


CO 


6^ Y 9u cl 

a T 8Y & 2 



(6.13) 


Is the vortlclty. 

The crucial step in the analysis Is to choose the relation between the 
amplitude scale e and the wavelength scale o so that the nonlinear terms 
produce a critical layer velocity jump of the same order as the velocity jump 
due to linear effects, l.e., 0(e). A little experimentation shows that 

nonlinear effects will influence the 0(e) term In u (both through the cross 
stream derivative In D and the vorticity source term on the right side of 
(6.10)) if we take 
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E = 


(6.14) 


Then u 1 and x 1 will satisfy 


(— - = -a — + b( Y - D /?e1A + e iX 

V3Y Vl J L3Y J 


Q) i + b(y 


1 )Y ) = 0 


where we have put 


b?--4=- /We iX U- 

V c l J ax 


^ \ N all 




/We'* iA f 


r-W e'H3- 

c ] ba l / 8Y 


This can be greatly simplified and put Into a more standard form by taking x ] , 


X, and 


Y = -Y + 

o bT c «,c, 


&AV* 


(6.18) 


as new independent variables. In which case ® becomes 


< 3 ) » aT + “ bY o — - -p- f^iA + e ix ) 

3x i 0 ax b V J 3Y 0 


(6.19) 


It now follows from (3.5), (4.13), (4.14), and (6.9) that 


2ir „® 


0 0 -® 


e" 1X — dY fl dX = - r4- A + ' , 

3Yn ° C l*i r 


( 6 . 20 ) 
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where we have put 


1 

r 


1 

T c q 


a. 



( 6 . 21 ) 


To put these results In a more universal form, we introduce the new 


variables 


H E 


( Y - nr 


Re aV x 

1 T c a l c l 


Q E - 


Vi 

br 


3u, 


9Y, 


“l( T c'l) 2 


Re A*. 1 *! 


t ,x o 

T A T e u 

A E ~ — s , 


bY, 


Y E 


r • 


X E X - X 0 , 


and 

* = r Vl - *0' 

Then H and 0 satisfy the homogeneous boundary conditions 

H.Q - 0 as Y - ±®, 
and (6.15), (6.16), (6.19), and (6.20) become 

0Q.(l-r + r ff)/?e1Ae 1X , 

S>h - -/?elAe 1X , 


and 


p® p2ir 
—CD 0 



dY dX 


dA 

dx’ 


( 6 . 22 ) 

(6.23) 

(6.24) 

(6.25) 

(6.26) 

(6.27) 

(6.28) 

(6.29) 

(6.30) 

(6.31) 
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where we have put 


» s Jr Y fe-(fc ,A,,x )J? (6 - 32) 

and 

r = (y - 1 >C 1 = T c - 1. <6-33) 

It is easy to see from these equations that 

H - -Re^z, <6.34) 

Y - 1 K 

and 

a- (\ - r)Re-^— <6.35) 

Y - lie 

when 

A - e K * as x - -®, (6.36) 

where 

k = (1 - r)ir (6.37) 

Is the scaled (and normalized) linear growth rate. The solution to these 
equations can therefore be made to satisfy the upstream matching condition 
(2.48) if we choose the, as yet, unspecified real constants Xq and xq to be 

X Q = -arg a, (6.38) 


and 


x o-^ ln Tir 


(6.39) 


7. NUMERICAL COMPUTATIONS 

The coupled nonlinear evolution equations (6.29) and (6.30) must be solved 
numerically. It is easy to see from these equations and the upstream boundary 
condition (6.36) that the solution A will remain real for all values of x. 
Then since Q and H are periodic in X, we expand them in Fourier series 
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J t V 


1 nX 


(7.1) 


h=- 


H - I E V* 


nX 


(7.2) 


h=-® 


with Q = Q*. H * -H* (where the asterisk, denotes the complex conjugate) to 
-n n -n n 


obtain 


(g; * ,nY ) n n * 2 A ay K+! * rH ntl " fl n-l " rH n-l] * ’ 5 n,i n ' r)A ' 

(fe * ,nY X * I A If < H n,l * Vl> - -'V, A - 


where A satisfies (6.36), 


S n i A 
H ^ 


n Y - Ik’ 


(1 - r)S .A 

aJ_ 


as x - -®, and 


- _ 

n Y - Ik 


9m \ a, 


dA 

dx 


(7.3) 

(7.4) 


(7.5) 


(7.6) 


(7.7) 


We solved (7.3) to (7.7) numerically using a procedure similar to the one 

used by Hanes (1985) and Goldstein, Durbin, and Leib (1987). Rather than 

mapping the Infinite domain -® < Y < ® Into a finite domain, they simply 

solved the equations over a finite range, say -N < Y < N, and used the 
asymptotic behavior of the solutions at Y - ±® to obtain an accurate 
approximation to the cross stream Integral ((7.7) in the present case). 

By using (6.29) and (6.30) to generate asymptotic expansions It Is easy 
to show that 


Hq.Qq 


0(Y" 3 ) 


(7.8) 
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(7.9) 


(1 - r) . 1(1 - r) dA „,„-3x 

Q, - w A + 5 — + 0(Y ) 

1 Y y 2 dx 


h. - - i a - i, ^ ♦ o<r 3 > 

1 T dx 


ii 


for n > 2 


as y - ±®. 


It follows that (7.7) can be approximated by 

[’ - I* 1 - r> ]gt - A. j" N 0, dy.OdC 3 ) 


(7.10) 

(7.11) 


(7.12) 


The numerical technique used to solve (7.3), (7.4), and (7.12) was similar 
to that used by Goldstein and Hultgren (1988). Equations (7.3) and (7.4) were 

solved for 0 < n < h where n was chosen so that the maximum absolute 

values of and Qfj remained below a certain present tolerance. The Y 
derivatives in (7.3) and (7.4) were discretized using second order central 
difference approximations and Simpson's one-third rule was used for the 
integrals in (7.12). The calculation was started In the upstream linear region 
where A, £) n , and H n are accurately approximated by (7.5) and (7.6) and 
marched forward In x through a predictor corrector procedure. A third order 

scheme was used for (7.12). The predictor step for (7.3) and (7.4) consisted 

of an Adams-Moul ton-Bashforth second order explicit approximation for the 
nonlinear and Inhomogeneous terms and the Crank-NIchol son (second order 
Implicit) approximation for the mean flow convection terms. The corrector step 

i 

was fully of the Crank-NI chol son type. The combined corrector steps for 
(7.3), (7.4), and (7.12) were then iterated until the solution at the next 
streamwlse station was obtained to within a preset tolerance. 
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8. NUMERICAL RESULTS AND DISCUSSION 

Goldstein, Durbin, and Le 1 b (1987) considered the nonlinear evolution of a 
two-dimensional Instability wave In a weak adverse pressure gradient boundary 
layer. Equation (6.29) reduces to their result In the limit where r - 0 and 
equations (6.29) and (6.30) become decoupled. They found that nonlinear 
effects always reduce the growth rate of the linear Instability wave, driving 
It toward an equilibrium state. 

The scaled Instability wave amplitude Is plotted as a function of the 
scaled and normalized streamwise coordinate x for various values of r In 
figure 1. The corresponding Instability wave growth rates, A^/A, are shown In 
figure 2. As In Goldstein et al . (1987) the growth rates Initially follow the 
linear growth until the amplitude becomes large enough for nonlinear effects 
to come Into play, but now the nonlinear effects can cause the growth rate to 
become larger than the linear growth presumably because compressibility 
effects, l.e., the Bjerknes forces, act as a vortlclty source within the 
critical layer. Notice that the growth augmentation Increases with 
Increasing r and only occurs when r exceeds a certain finite value, say 
r o- 

Figure 3 shows the rollup of the vortlclty contours In the X-Y plane at 
various (Increasing) values of x. The rollup Is not too different from that 
found by Goldstein et al . with the principal difference being the formation of 
an additional counter rotating vortex core. The results are replotted In 
figure 4 versus the more physical coordinates X and -Yb/r for 

r a i 
Wi 2 

Notice that they now look considerably different. 
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